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■ In this paper the SU{N) Einstein-Skyrme system is considered. We express the chiral field (which 



is not a simple embedding of the SU (2) one) in terms of harmonic maps. In this way, SU{N) 
spherical symmetric equations can be obtained easily for any N and the gravitating skyrmion 
solutions of these equations can be studied. In particular, the SU (3) case is considered in detail 
and three different types of gravitating skyrmions with topological charge 4, 2 and 0, respectively, 
are constructed numerically. Note that the configurations with topological charge correspond to 
mixtures of skyrmions and antiskyrmions. 
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I. INTRODUCTION 



Nonlinear field theories coupled to gravity have received a lot of interest in the past 
decade. It has been discovered that gravitational interaction may lead to genuinely nonper- 
turbative phenomena like gravitationally bound configurations of nonabelian fields, while 
the study of black hole solutions in various models revealed the possibility of nonlinear hair 
on black holes which questioned the validity of the unqualified no-hair conjecture. 

One of the candidates of such investigations is the Einstein-Skyrme model which can 
e.g. describe the interaction between a baryon and a black hole (a configuration which 
might have been produced in the very early universe). So far, most of the studies have 
concentrated on the SU{2) Einstein-Skyrme model Ql- In particular, in it has 

been shown that the Schwarzschild black hole can support chiral ("Skyrme") hair and it 
has been argued that such configurations might be stable. The presence of the horizon in 
the core of the skyrmion unwinds the skyrmion, leaving fractional baryon charge outside 
the horizon. More systematic investigations of the model were undertaken in j^l by solving 
numerically the static spherically symmetric Einstein-Skyrme equations. Globally regular 



solutions with baryon number one P] and black holes with chiral hair were found and their 
stability properties were studied in detail 0,0]. 

The first examples of nonembedded solutions for a higher group, namely the SU{3) 
group, were the S0{3) solitons with even topological charge. The lowest energy solution 
corresponds to a bound state of two gravitating skyrmions Specifically, it was found 
that the two branches of these regular solutions exist, which merge at a critical value of the 
gravitational coupling. 

In this paper we consider particle-like solutions of the SU{N) Einstein- Skyrme model 
for N > 2). In particular, we study the deformation of the multiskyrmion configurations 
^ (derived using the harmonic map ansatz) when gravity is introduced. Note that the use 
of harmonic maps for gravitating skyrmions can be traced to Ql- New types of solutions 
are found, which correspond to skyrmion-antiskyrmion configurations and have topological 
charge 0. Like the non-gravitating skyrmion-antiskyrmions, these configurations are also 
saddle points of the energy functional and thus are likely to be unstable. 

Our paper is organised as follows: in Section II, we present the SU{N) Einstein-Skyrme 
model, while in Section III we give the harmonic map ansatz. In Section IV, we present the 
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spherically symmetric equations of motion and in Section V we discuss our numerical results 
for the SU (3) case. In this latter section, we also point out that we recover the already known 
solutions for the topological charge 4 and 2 jsl, but also derive new solutions describing 
the gravitating skyrmion-antiskyrmions with the topological charge 0. Our conclusions are 
summarised in Section VI. 



II. THE SU{N) EINSTEIN-SKYRME MODEL 

The SU{N) Einstein-Skyrme action reads: 
R 



S 



tr (K^ K^") - — tr i\K^, KJ IK^", K""]) 



-g d'^x 



(1) 



where = d^UU^^ for /i = 0, 1, 2, 3 and f/(x^) G SU{N) is the matter field, g denotes the 
determinant of the metric and G represents Newton's constant. In order for the finite-energy 
configurations to exist the Skyrme field has to go to a constant matrix at spatial infinity: 



U —>■ I a.s \x^ 



oo. 



To derive the classical equations of motion of our system we perform the variation of the 
action (0) with respect to the metric and the Skyrme field. The variation with respect to 
the metric g'^'^ gives the Einstein equations: 



(2) 



where R^^ denotes the Ricci tensor and the stress-energy tensor T^^, = g^u^ — 2^^ is given 
by: 

T^,= tr (^K^K, - ]^g^,K^K^^ + ^tr (^g^" [K^, K^] [K,, Kp] - ^g,, [K^, K^] [iT", K"] ^ . 

(3) 

The variation of the action with respect to the matter fields leads to the Euler-Lagrange 
equations which we will discuss in the next section. 

The Einstein-Skyrme systeni__has a topological current which is covariantly conserved, 
yielding the topological charge 



B 



-g d^x 



(4) 



where 



B^ 



247r2, 



(5) 



and e'^^"^ is the (constant) fully antisymmetric tensor. 

In what follows we will concentrate our attention on studying the static Einstein- Skyrme 
equations and we are going to construct their static spherically symmetric solutions by using 
harmonic maps for the Skyrme field. 



III. HARMONIC MAP ANSATZ 



The starting point of our investigation is the introduction of the coordinates r, z, z on 



K^. In terms of the usual spherical coordinates r, ^, the Riemann sphere variable z is given 
by: z = e*'^tan(^/2) and z is the complex conjugate of z. In this system of coordinates the 
Schwarzschild-like (spherically symmetric) metric reads: 



1 



Cir) 



'l + \z 



2\2 



dzdz, C*(r) 



1 - 



2m{r) 



(6) 



where m(r) is the mass function. For this metric the square-root of the determinant takes 
the simple form: 



2r' 



(1+|Z|2)2 

After substituting the metric ©, the action (0) becomes: 



(7) 



S 



R 



16nG 



tr 



1(1 + 1 



\2\2 



2\4 



2r2 



\2\2 



-c 



2 



1 (1+ 1^ 

32 

dt dr dzdz 



while the baryon number is equal to 



tr {Kr [Kz, K2]) dr dz dz 



In addition, the Einstein equations take the form: 



—m 



87rGT° 



2 A' 
r A 



c = sttG (T° - t;) 



where the prime denotes the derivative with respect to r and 



J n 



1 (1+ 1^ 

'32 r4 



2r2 



tr 



2\4 



T° - T; = -Ctr {K^) - l^^i^J^tr 



(8) 



(9) 



(10) 



fill 



(12) 



Following |6[ the application of the harmonic map ansatz to describe the matter fields, 
corresponds to setting: 



= e-^i9o/N ^ ^-2ig,/N + A^P,)... (1 + P^_2) (13) 

where (7^ = (7a: (r) for k = 0, . . . , N — 2 are the profile functions which depend only on r. 
Moreover, we define also Ak = e^*^*-' — 1. The boundary value U ^ / at r ^ 00 (needed for 
finiteness of the action) imposes the requirement that gi{oo) = 0. 

Here Pk form a set of projectors based on the maps S"^ CP^^^ which are constructed 
as follows 0]: write each projector P as 

V 

nv) = (14) 

where \^ is a A^-component complex vector of two variables z and z. The first projector is 
obtained by taking V = f{z) (i.e. an analytic vector of z), while the other projectors are 
given in terms of new vectors V which are obtained from the original V by differentiation 
and Gramm-Schmidt orthogonalisation. If we define an operator P+ by its action on any 
as 



vector V E (E 



N 



3 



P+v = d,v-v I , (15) 
\v\ 



then the vectors = P+v can be defined by induction: P^v = P+(Pf;^^f). 

Therefore, in general, we can consider projectors Pk of the form ()14|) corresponding to 
the family of vectors V = Vk = P^f (for f = f{z)) a.s 

Pk = PiPlf), k = 0,...,N-l, (16) 

where, due to the orthogonality of the projectors, we have "^^=0 = 1- This follows from 
the following properties of vectors P^f (which hold when / is holomorphic) : 



g_(pkf) ^ _pk-ifjf_±lj_ g _Lj_J_] = _L±J_ (ij) 



(P^/)tp|/ = for M/, 

\plf? . ( pT'f \ ^ Plf 

Note that, for SU{N), the last projector in the sequence, Pn-i, corresponds to an anti- 
analytic vector (i.e. a function of z). Moreover, we can always express one projector as a 
sum of the others. 
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It was shown in p] (for the non-gravitating spherical symmetric skyrmions) that the 
chiral field (|T^ is an exact solution of the corresponding equations when 

f={fo,...,f„...,fN-iy where f^ = z^^(^^J^^ (18) 

and {^J^) denote the binomial coefficients. In what follows we apply the same ansatz to 
obtain the corresponding gravitating skyrmions for the simplest cases of SU{2) and 5*^/(3). 
These cases will clarify the expressions for the general SU (N) case. 



IV. SPHERICAL SYMMETRIC EQUATIONS OF MOTION 



In the case of spherical symmetry, the action ([Q) using (jl3j) takes the form 



S = 271 



+ 



A 

47-2 



RAr^ 4 
8nG ~ N 

N-2 



m-2 



N-2 



N-1 



'CAiJ^g'A +4r'CAj29? + ^Aj2D^ 



. i=0 



1=1 



i=0 
N-l 



k=l 



+ 



2CAY,D,{9',-9',_,y]drdt, (19) 



k=l 



where = 2k{N — k) siia^{gk — Qk-i)- The matter equations are obtained from the vari- 
ation of this action with respect to the matter field. We will present these equations when 
considering the specific cases of SU{2) and SU{?>). 

In addition, the Einstein equations (fTUI) take the form: 



'N-2 



N-2 



N-l 



N-l 



—m' = 16ttG 



. i=0 
N-2 



1=0 



2r' 



k=l 



2r' 



k=l 



16r^ 



+ (A - A+i)' + Dl_, 



2 A' 
r A 



C = 16nG 



2c; 



i=l 



'N-2 



N-2 



N-l 



J2g'A +2Cj2g? + C^^J2D^U-9Li) 



i=0 



i=0 



k=l 



(20) 



(21) 



For simplicity, we set Fj. = — g^^i for k = 0, . . . , N — 2 with -F/v-2 = gN-2- Then, the 
topological charge ^ simplifies to 

B='-J:i^ + l)iN-^-l)(F^-'^) . (22) 



Since Fi{oo) = the only contributions to the topological charge come from -Fj(O). 
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The main symmetry of our expressions is the symmetry under the independent inter- 
changes 

k = 0,...,N-2 (23) 

which follows from the fact that the terms in ()19p are symmetric under the interchange 
Dk ^ D^-k when Fk-i ^ F]\f^k-i- At the same time, all the other terms in exhibit 
this symmetry since they are combinations of Fi and their derivatives. 



A. SU{2) 



For N = 2 there is only one profile function, Fo{r), and (fT^ simplifies to 
S = 2ti I [ + 2r^CAF'^ + AA sin^ Fq + 2A^^^ + AC A sin^ FqF'^ ) drdt 



8nG 



while its variation with respect to the profile function gives us the matter equation: 



CAr^F(^ ( 1 + 



2 sin^ F, 



Asin2Fo 1 + 



0. 



The Einstein equations (j?H) now take the form: 
2 , 



SttG 



sin^ Fn 



sin^ Fo 



-4^ = 167rGF^2( 1 + 2C „ 
r A \ 



(24) 

(25) 

(26) 
(27) 



There is only one solution with the boundary condition FnfO) = tt. This solution has 

nn 

topological charge B = \ (see ((221)) and has previously been studied in great detail in j3|,y|. 
Therefore we will not repeat these calculations here. 



B. 5?7(3) 



For = 3 there are two profile functions, Fo(r), Fi(r), and (|19|) becomes 

^ = 2^ / I + ^ ^'^^ (^o' + ^i" + ^oi^O + 8^ (si^' ^0 + sin^ F^) 

+ ^ (sin^ Fo - sin^ Fq sin^ Fi + sin^ Fi) + %C A (sin^ FgFfJ^ + sin^ FiFf ) jdrdt (28) 



The corresponding equations for Fq and F\ are now given by: 



r^CAF' ; 



2 2sin2Fo\ 



-Asin2Fo 1 + 



2 sin^ Fo - sin^ Fi 



+ CF 



72 



(29) 



r'^CAF[ 



2 2 sin^ Fi 
7: + ^ 



-Asm2Fi 1 + 



2 sin^ Fi - sin^ Fq 



72 



(30) 



Note that the above equations are symmetric under the simultaneous interchange Fq —>■ Fi 
and Fi ^ Fq. 

Finally, the Einstein equations take the form: 

I (i^o' + + ) + ^ (sin^ Fo + sin^ F,) + ^ (sin^ FoF^^ + sin^ F^F^^) 



— m' = 327r G 



L3 



1 



-4^ = 647rG 



+ — (sin^ Fo - sin^ Fq sin^ Fi + sin^ Fi) 
1 (F- + F^F; + Ff ) + ^F- + ^Ff ■ 



(31) 
(32) 



The set of equations ()29p - (p?^ can only be solved numerically when the right boundary 
conditions have been imposed. Similarly to the flat case we see that there exist three 
types of gravitating multiskyrmions, which we will discuss in detail in the following section. 



V. NUMERICAL SIMULATIONS 

To solve the equations ()29|l - (jH^ numerically, we have adopted the numerical routine 
described in For convenience, we define = IGvrG so that the fiat limit with C{r) = 
A{r) = 1 corresponds to a = 0. 

A. Boundary conditions 

As will be discussed in more detail in the following, three kinds of regular finite-energy 
solutions exist with the following boundary conditions: 



(I) Fo(0) = TT, 

(II) Fo(0) = TT, 

(III) Fo(0) = TT, 



Fi(0) = vr, 

Fi(0) = -TT, 

Fi(0) = 0, 



(33) 



and for all these cases 



(I) / (II) / (III) : Fo(oo) = , Fi(oo) = 0. (34) 
Furthermore, we have the two supplementary conditions for the metric functions: 

(I) / (II) / (III) : m(0) = , A{oo) = 1 . (35) 
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The condition that m(r) vanishes at the origin r = is due to the requirement of regularity, 
while the second condition for A{r) results from the requirement of asymptotic flatness. The 
energy E of the gravitating skyrmions can then be determined from the "mass function" 
m(r) at infinity: 



E = . (36) 



4m(oo) 

With this normalisation, the values of E can be compared to those of the flat limit 

Our numerical analysis demonstrates that the three solutions are indeed continuously 
deformed by gravity (a > 0). Next we discuss our numerical results for each set of conditions 
given in 

B. Case I 

This case corresponds to choosing -Fo(r) = Fi{r). Apart from a trivial rescaling of the 
coupling constants p|, the equations for the matter and metric functions are equal to those 



a 



However, 



of the gravitating SU{2) skyrmion, which were studied in great detail in 
for completeness, we again present the main features of these solutions. 

The non-gravitating solution has energy E ^ 4.928 = 4 x 1.232, i.e. four times the 
energy of the SU{2) one-skyrmion. Due to the boundary conditions -Fo(O) = -^i(O) = vr the 
topological charge is four (see (j221))- The non-gravitating solution can thus be interpreted 
as four noninteracting skyrmions placed on top of each other in such a way that the baryon 
(energy) density is spherically symmetric. 

As for the asymptotic behaviour, we notice that the chiral fleld is of the form 

Fo{r) = Fi{r) ^ n - Bjr for < r << 1 , 

F^{r) = Fi{r) ^ for r >> 1 , (37) 

where 5/, Bi are (shooting) parameters depending on a which have to be determined 
numerically. 

Solving numerically the system ()29|l - (jH^ for a > we flnd that the flat solution is 
gradually deformed by gravity, forming a branch of gravitating skyrmions. In particular, the 
function C(r) develops a local minimum at some intermediate radius : r = Vmia), while the 
function A{r) has a minimum A^in = A{0) at the origin and then increases monotonically. 
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However, the profile functions of the Skyrme field deviate only slightly from the corre- 
sponding ones in the fiat limit. Moreover, as a increases, the respective minimal values of 
the metric functions C(r), A{r\ i.e. Cm = C{rm) and Aq = A{0) both decrease and so does 
the corresponding energy E. This latter decrease is, of course, expected since gravity tends 
to lower the mass of a solution. The four skyrmions can thus be seen as "gravitationally 
bound". This is illustrated in Fig. ^ (solid lines). Note that this branch of gravitating 
skyrmions does not exist for arbitrarily large values of a, but only up to some critical value 
o^cr- Oi < acr ~ 0.142087 (as shown in Fig. P). Also the quantities Aq, Cm as functions of 
a remain finite with 

E{a = a^r) ^ 4.20, Ao(a = a^r) ~ 0.437, Cm{a = a^r) ^ 0.584. (38) 

Our numerical analysis, however, strongly suggests that a second branch of solutions 
exists in the interval [0,acr]- For a given a, the solution of the second branch has a higher 
energy than the one on the main branch while Aq, Cm have lower values as shown in Fig. 
H (dotted lines). For a acr both branches go to the same solution but the solution of 
the second branch becomes more and more peaked around the origin (i.e. the slope F'{0) 
tends to infinity). As a consequence (in this limit) the energy on the upper branch diverges 
as a — ^ 0, but the product aE remains finite (~ 0.124). This solution, rescaled according 
to X ^ - stays regular in the a — >■ limit and converges to the sphaleron solution of 

nn 

the Einstein- Yang-Mills system |5|, 112|. Note that the metric functions remain finite and 
are strictly positive. Therefore, no black hole solution is generated by the solutions of the 
equations under consideration, in contrast to e.g. the Einstein- Yang-Mills-Higgs equations 



C. Case II 



By choosing Fo{r) = —Fi{r) with the boundary conditions -Fo(O) = — -Fi(O) = vr, the 
nongravitating solution is topologically trivial since the topological charge is zero. How- 
ever, this configuration is not the vacuum solution, but consists of two skyrmions and two 
antiskyrmions. Note that these solutions were not studied previously in jSj]. 

The regularity of the solutions at the origin requires that Fq{0) = which has been 
checked by us numerically to hold within a very high level of accuracy. The flat solution 
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(a = 0) has energy E ^ 3.861 and was discussed in detail in p. 

It is worth noticing that the asymptotic behaviour of the Skyrme function -Fo(r) differs 
drastically from the one of case I: 



Fo(r) ^ vr - -B//r^ for < r << 1 , 

Fo(r) ^ for r»l . (39) 

For a > the pattern of the solutions is very similar to the one occuring in case I. The 
critical value of a is larger than in I with a^r ~ 0.1834. This is related to the fact that the 
solutions of case I are heavier (compare the flat limit energies) and thus exist for a smaller 
interval of the gravitational coupling. We note that 

E{a = «„•) ~ 3.377, Ao{a = a^r) ~ 0.517, Cm{a = a^r) ~ 0.614. (40) 

The occurence of two branches is illustrated in Fig. 2, where the energy E of the solution 
and the quantities Cm, are plotted as a function of a. 

In the limit a ^ on the upper branch, the solution again converge to the SU(3) 
Einstein- Yang-Mills sphaleron solution constructed first in . When solving the equations 
numerically, the main feature of the present charge solution resides in the fact that the 
derivative of the chiral function -Fo(r), E^lr), vanishes at the origin r = and thus the 
"shooting" parameter is the second derivative FQ{r)\r=o- This renders the numerical analysis 
more involved than in the first case. The solutions' profiles corresponding to the lower and 
upper branches are shown in Fig. 3 for a = 0.17. In this figure, a logarithmic scale is 
used in order to reveal the completely different behaviour of -Fo(r) as r ^ 0. The value 
FQ{r)\r=o for the solution on the upper branch is several orders of magnitude larger than its 
counterpart on the lower branch. More generally, we noticed that the shooting parameter 
FQ{r)\r=Q varies surprisingly strongly with a and increases considerably when a increases 
(resp. decreases) on the lower (resp. upper) branch. 

D. Case III 

In this case we choose Fq{0) = tt, -Fi(O) = and Fi{r) goes from zero to zero developping 
one node at some finite value of r. Due to the boundary conditions, these solutions have 
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topological charge 2 and are the gravitating 5*0(3) embedded solutions discussed in [5|. The 
flat limit {a = 0) solutions have been studied in Q]. Again, for completeness, we discuss 
the main features of the solutions in the following. 

The regularity of the solutions at the origin implies that Fq{0) = F[{0), which confirms 
the validity of our numerical procedures. The energy of the fiat solution is ~ 2.376. The 
pattern with two branches of solutions merging at a critical value of acr ~ 0.2333 is also 
seen here j^. 

The gravitating solution can be characterized by 

E{a = acr) ~ 2.033, Ao{a = acr) ~ 0.460, Cm{a = acr) ~ 0.537. (41) 



VI. CONCLUSIONS 



Our numerical results (for the three cases studied) indicate that the higher the energy 
of the solution is, in the fiat limit, the smaller is the interval of a for which the gravitating 
skyrmions exist. The metric functions of the solutions on the second branch do not develop a 
zero in the critical limit and thus do not represent solutions with horizons. It is rather that, 
in this limit, the matter fields become singular at the origin. However, using a different scale 
for the radial variable (e.g. by reinstating the Skyrme coupling constant) would render the 
limiting solutions regular. Following jl^ , we call such solutions "gravitating sphalerons" . 

Following the arguments based on Morse theory (see eg and references therein), we 
note that the solutions on the branch with the higher energy have one more unstable mode 
than the solutions on the lower branch. In |5] it was argued that the solutions on the lower 
branches in cases I and III are stable (since the fiat space limit on these branches is stable), 
while the solutions on the upper branch are unstable. For case II, the fiat space solution 
on the lower branch is unstable and leads to annihilation P|. Thus, we expect that the 
gravitating analogues on the lower branch are unstable and that the solutions on the upper 
branch have two unstable modes. 

Let us mention that the situation discovered here for the cases I and III is, in many 



m, Ua. In the latter 



respects, similar to what happens with the electroweak skyrmion 
case the starting theory is the standard model of electroweak interactions, in which the 
Higgs field plays the role of a Skyrme field and the relevant coupling constant (call it ^) 
parametrizes a supplementary effective interaction (a sort of Skyrme term) encoding part of 
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the radiative corrections of the theory. Two branches of solutions exist which also merge at 
a critical value of C,- The lower branch is stable and is the so-called "electroweak skyrmion 
branch" because it approaches a skyrmion solution in the limit ^ — » 0. The higher branch, 



as shown in 17|, approaches either the sphaleron 
Higgs field mass) the bisphaleron solution jiol . 2^ . 



or (according to the value of the 
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FIG. 1: The values of the metric function A(r) at the origin A{0) = Aq, the minimal value of 
the metric function C (r) , Cm and the energy E = ^ "^j^^ of the gravitating skyrmion solutions 
for case I are plotted as functions of a. The solid (respectively dotted) lines refer to the first 
(respectively second) branch of solutions. 
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FIG. 2: The values of the metric function A(r) at the origin A{0) = Aq, the minimal value of 
the metric function C (r) , Cm and the energy E = ^ "^j^^ of the gravitating skyrmion solutions 
for case II are plotted as functions of a. The solid (respectively dotted) lines refer to the first 
(respectively second) branch of solutions. 
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FIG. 3: The profiles of the function Fo(r), FQ{r) and C(r) are shown for the two branches of 
skyrmion-antiskyrmion solutions (case II) with a = 0.17. 
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